Dipole perturbations of the Reissner-Nordstrom solution: II. The axial case 
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We study the linear metric perturbations of the Reissner-Nordstrom solution for the case of axial 
perturbation modes. We find that the well-known perturbative analysis fails for the case of dipole 
(I = 1) perturbations, although valid for higher multipoles. We define new radial functions, with 
which the perturbation formalism is generalized to all multipole orders, including the dipole. We 
then complete the solution by constructing the perturbed metric and Maxwell tensors. 
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I. INTRODUCTION 

The linear gravitational and electromagnetic pertur- 
bations of the Reissner-Nordstrom solution have been 
analyzed by Moncrief ||l|, Zerilh [Q, Chandrasekhar Q, 
and Chandrasekhar and Xanthopoulos Q. (For a sum- 
mary of these works see Refs. [||^. These perturbations 
were considered also in 0-0 •) In the perturbative ap- 
proach the deviations from the exact Reissner-Nordstrom 
solution are considered as (small) perturbations over the 
Reissner-Nordstrom background. Due to the coupling of 
the electromagnetic and the gravitational fields in the 
background, any gravitational perturbation leads in gen- 
eral via the field equations to an electromagnetic per- 
turbation, and vice versa. Despite this complication it 
turns out that the equations describing the evolution of 
the perturbations can be decoupled (for each multipole 
order and for each parity) to a pair of independent one- 
dimensional wave equations (namely, the Moncrief- Zerilli 
equations), each of which is made of an electromagnetic 
component and a gravitational component. 

In a recent paper |ll| (hereafter Paper I) we analyzed 
the polar modes of the perturbations, and showed that 
the perturbative formalism of Ref. [^ fails in the case of 
dipole perturbations (namely, the I — 1 modes in a multi- 
pole expansion). We generalized the perturbative analy- 
sis to include also the dipole modes by defining new radial 
functions, and completed the solution by giving the met- 
ric perturbations and the perturbations of the Maxwell 
tensor in terms of the perturbing fields. In this sequel we 
treat the axial perturbations of the Reissner-Nordstrom 
solution. (For a definition of the parities see, e.g., Ref. 
P|.) As we shall see in what follows, also in the axial case 
the formalism of Ref. H fails, and for similar reasons to 
the failure in the polar case discussed in Ref. h% . In this 
paper we shall again define new radial functions, with 
which one can formulate a perturbative analysis which is 
valid for all multipole orders, including the dipole. 

For many applications the failure of the perturbative 
formalism for the dipole modes is immaterial. This is 
due to the non-radiative character of dipole gravitational 
modes. However, there are cases where dipole electro- 
magnetic waves are important, especially in the Reissner- 



Nordstrom spacetime, because of the coupling of the elec- 
tromagnetic and gravitational fields. For example, the 
late-time behavior of the electromagnetic perturbations 
produced during the collapse is controlled by the dipole 
modes (which are the least-damped modes) |ia-H]. In 
addition, the dipole perturbations are especially impor- 
tant in the analysis of the (electromagnetic) perturba- 
tions of the Cauchy horizon of Reissner-Nordstrom black 
holes |13]. One needs, therefore, a valid perturbative 
formalism for the calculation of the perturbed metric co- 
efficients and the components of the Maxwell tensor. 

The organization of this paper is as follows: In sec- 
tion II we discuss our definitions and the notations. In 
section III we review the perturbative approach we use. 
In section IV we derive the coupled equations, which we 
decouple in section V, and in section VI we complete 
the solution, which is adequate for the treatment of the 
dipole modes. 



II. DEFINITIONS AND NOTATION 

Following the notation and convention of Ref. [||, 
we write the line element of an unperturbed Reissner- 
Nordstrom black hole in the form 



ds2 = e"" ( da:°) ' - e^^^ ( dx'f ~ r^ dVi} 



(1) 



where d^ 



o2V' 



dx^y 



s^'^3 {dx^\ is the line el- 



ement on the unit two-sphere, the metric coefficients 



o2f 



-2^2 



(^2 _ 2Mr -f Ql)lr^ 



^jr^ 



and e^"^ = e^^^ sin^ 9 — r^ sin^ 9, and the coordinates 
are {xP x^ x^ x"^) = {t (f> r 9) . Here, r is the radial 
Schwarzschild coordinate, defined such that circles of ra- 
dius r have circumference 27rr, and M and Q* are the 
mass and electric charge, correspondingly, of the black 
hole. 

It turns out that the metric perturbations of the 
Reissner-Nordstrom black hole can be separated to two 
parities — called polar (or even parity) and axial (or odd 
parity) — in accordance with their behavior under the 
transformation (f> ^ —(f>. The Reissner-Nordstrom back- 
ground, being spherically symmetric and static, is charac- 
terized by the vanishing of all the axial metric coefficients 



[namely, in the background one has w = (72 = 93 = 0; see 
the definitions for these metric coefficients below in Eq. 
(0)]. Thus, any axial perturbation will impart rotation 
on the black hole, and transform it into a (slowly) rotat- 
ing Kerr-Newman black hole. It can be shown H that 
the metric 



ds^ = e^-" {dx°Y - e"^ {dx' - u dx"" - 92 dx' 
- q^dx'f - e^f^- {dx^f - e^^^' {dx'f . 



(2) 



is of sufficient generality for the treatment of all pertur- 
bations (polar and axial). 

Throughout this paper we shall use the notation of Ref. 
H, except when we change its formalism; Then, simi- 
larly to the notation of Paper I, we shall add a bar to 
the symbols of Ref. Q. The 'barred' objects will be de- 
fined in such a way that dipole perturbations are treated 
adequately. This notation facilitates the comparison be- 
tween our formalism, and the formalism of Ref. IHI. 



III. THE PERTURBATIVE FORMALISM 

In a similar way to the development of the perturba- 
tive equations for the polar modes, the equations gov- 
erning the axial perturbations are derived by a lineariza- 
tion of the coupled Einstein-Maxwell equations about the 
Reissner-Nordstrom background 0. In particular, after 
substitution of the expressions for the unperturbed met- 
ric coefficients (of the background) in the linearized (j) f 
and (f> 9 components of the Ricci tensor, one obtains the 
equations 

[r^e^''{q2,3 - 93,2) sin^ 6] 3 + r^{uj.2 - 92.0), sin^ Q 



= 2r3e''sin2 6'5i?(i)(2) 

= 4g,re''F(o)(i) sin^ Q (3) 



and 



{r^(?\<i2,z - 93,2) sin^ e\ 2 - r2e-2"^(a;,3 - 93,o),o sin^ B 

= -2r2sin2 6'(5i?(i)(3) 
= 0, (4) 

where Xi^\{0\ is the (a)(/3) tetrad component of the ten- 
sor X. Here, R and F are the Ricci and Maxwell tensors, 
respectively, and a comma denotes partial differentiation. 
Substitution in the linearized components of the Maxwell 
tensor yields 



h"^(0)(i)],2 + ^e "^(i)(2),o = 0, 



(5) 



[F(o)(i) sin B\ 3 -f r2F(i)(3)^o sin0 = 0, (6) 



re "^(o)(i),o+ [?'e"^(i)(2)]_2 + ^(i)(3),3 

= -(5*(tJ2 - 172,0) sin6'. (7) 



We now introduce the functions 



B(r, 



Q{r,9) = 



F, 



(0)(1) 



'(92,3 



'73,2) sm 



(8) 



(9) 



After substitution of these functions in Eqs. (||) and fl), 
and using the linearized Maxwell equations [Eqs. (0)- 
(0)], we obtain the differential equations 



dQ 



r4o;r,3 5» d0 



r* sm 



-(^,2 -g2,o).o + 



4g* 



-e^B 



and 



dQ 



r^ sin^ ( 



dr 



{^,: 



93,0J,0- 



(10) 



(11) 



We now assume that the perturbations can be analyzed 
to their normal modes with a time dependence e*°'* . (This 
is always possible for a linear theory.) Because of this 
explicit time-dependence of the metric perturbations we 
did not consider the time dependence in the definitions 
for the functions B{r, 9) and (5(r, 9). Substitution in Eqs. 
(|l^) and (|ll|) yields then 



1 



dQ 



.4,;„3g, QQ 



r^ sm 



= —lauj 2 



cr^q2 



4Q* 



r^ sin'^ 



-e^B (12) 



and 



A 



dQ 



1^9 dr 



iauj^s + cr (73. 



(13) 



Differentiating Eq. (h2[) with respect to 9 and Eq. (nj 
with respect to r, and summing the equations, we find 
that 



d /A dQ 



dr V r^ dr 



sm 



_d_ 
'd9 



1 dQ 
n^9 d9 



= 4Q.e^.A 



B 



sm 



sm 



(14) 



We now differentiate Eq. (|) with respect to r, Eq. (|g) 
with respect to 9, and Eq. (TJ) with respect to t. Substi- 
tuting in the latter and using Eq. m), we find 



[e2-(re''i3),2] , + 



B' 



smf 



sinf? — re '^ B qq 



Q*(w,2,o - 92,0,0) sin^ 



(15) 



Eqs. (|14| ) and (|15[) govern the axial perturbations. 

We next separate the variables in B{r, 9) and Q{r, 
This separation is done in Ref. |5[| by the ansatz 



and 



B{r,9)=iB{r)C 1^1^ {cob 9) 



Q{r, 9) ^Q{r)C;_^l^ [cos 9) 



(16) 
(17) 



where C^ is the Gegenbauer function of order n and in- 
dex V |1^. It turns out that the relations given in Ref. 
[gl between the Gegenbauer functions and the Legendre 
functions are correct only for I > 2. (See Eqs. (21) 
and (22) in Chapter 4 of Ref. g.) For I = 1, how- 
ever, they give erroneous results. The correct expression 
for the n = 3 and v — —3/2 Gegenbauer function is 
Cg" (cos 0) = ^ (cos^ — 3 cos 9) , whereas the expres- 
sions given in Ref. S give (incorrectly) an identically- 
vanishing expression. For reference, we give here an ex- 
pression from which the Gegenbauer functions can be 
calculated directly for all orders and indices fig] : 

r^^f a\ v^ r(z/ -h m)r(j/ + n - m) 

C;^(cos6') = > — cos{n~2m)0. 

^0 m!(n-m)![r(z.)]2 

Substituting Eqs. (|l|) and (0) in Eqs. (Q and (|l|) 
we obtain the radial equations 



dr \r^ dr I ^^ 



A 



^_1^(/ _!)(; + 2) Ae'-B 



(18) 



and 



d 
dr 



e^'' — (re''B) 



(J re' 



dr 



l(^l + l)_B 



eM B = -g. 



Q 



(19) 



IV. DERIVATION OF THE COUPLED 
EQUATIONS 

We now define the functions 



H 



(-) 



and 



H. 



(-) 



-2re''B{r) 



-Q{r) 



(20) 



(21) 



divisions and multiplications of finite physical fields by 
identically-vanishing expressions, and thus obtain non- 
sensical results. It should be stressed, however, that for 
all other modes, namely, for / > 2, the perturbative for- 
malism as presented in Ref. 0] is perfectly correct and 
valid. 

We now change variables to the Regge- Wheeler "tor- 
toise" coordinate r* defined by dj dr^ = [I^/r'^)d/ dr, 



r(-) 



(-) 



and find that H{ and H^ satisfy a pair of coupled 
second-order differential equations 



A^ijf-) 



A 



/(/ + l)r-3M + 4 



Ql 



H 



(-) 



r(-) 



(-) 



and 



A^ijf-) = 



A 



3MHI ' + 2Q^H. 



l{l + l)r - MI + A- 



(22) 



H. 



(-) 



ZMH. 



(-) 



+ 2Q^{1-1){1 + 2)h[-'^Y (23) 



Here, A^ = d^ / dr^ + ct^. It is important to notice that 
Eqs. ( p2| ) and ( p3| ) are already decoupled in the following 
two cases: First, when Q, = 0, namely, when the elec- 
tric charge of the black hole vanishes and the Reissner- 
Nordstrom black hole degenerates to Schwarzschild. In 
this case there should be just one radiative perturbative 
mode, described by the Regge- Wheeler equation, which 
describes the perturbations of Schwarzschild. Indeed, we 
find that Eq. ( P3|) reduces to the Regge- Wheeler equa- 
tion in the limit of vanishing electric charge. The other 
equation apparently describes another radiative mode (as 
it is a wave-like differential equation). We shall discuss 
the meaning of the other apparent mode in the next sec- 
tion. Second, when 1 = 1, namely, in the case of dipole 
perturbations. (In fact, in this case Eq. (E3|) governing 



r(-) 



r(-) 



H2 is decoupled, while Eq. (g2|) still couples H^ and 
H2 ■ However, in such a case it is very easy to decouple 
the equations.) In the next section, after we decouple the 
equations, we shall discuss this case in detail. 



Note, that here we deviated from the formalism of Ref. 
M. The reason for this deviation is as follows: In Eq. 
(142) of Chapter 5 of Ref. 0] the variables /i and n (do 
not confuse n here with the order of the Gegenbauer func- 
tions) are defined by fi'^ = 2n = {l — l){l + 2). Then, when 
the function H\~' is subsequently defined (Eq. (143) of 
Chapter 5 of Ref. |g]) it is divided by /i. However, for 
dipole modes / = 1, and consequently fi vanishes iden- 
tically. As B and H[ are physically-meaningful func- 
tions [they determine the electromagnetic field via Eq. 
(0) and the Maxwell equations] this is manifestly inap- 
propriate. By keeping this ill-defined functions for the 
treatment of dipole perturbations, one might encounter 



V. THE DECOUPLING OF THE EQUATIONS 

A. Decoupling the equations for the dipole case 

In a similar way to the decoupling procedure we used 
in the polar case in Paper I, we first examine the / = 1 
case. As already noted in the preceding section, in this 
case the equation governing the evolution of H2 is al- 
ready decoupled from 7J| . It is easy to verify [by direct 
substitution in Eqs. ( |2^ ) and (p3|)] that the complete de- 
coupling can be obtained by a transformation to the new 



-A~) 



(-) 



radial functions Zl and Zj defined by 



H 



(-) 



91 



Z 



(-) 



.0. 

91 



K-) 



and 



iJ. 



(-) 



6M 
91 



(-) 



(24) 



(25) 



Here, gi is an arbitrary parameter. (It will be fixed when 
we decouple the equations for any I below.) The factor 
(— 6M) in Eq. (£3) was taken for convenience. We then 
find the equations for the functions Z\~' and Z^ to be 



A^Z^-^ 



A 



2r 



,Qt 



Z- 



{-) 



and 



A^Z, 



(-) 



= — ( 2r - 6M 



,91 



-A-) 



(26) 



(27) 



Eq. ( p7| ) is just the Regge- Wheeler equation for the case 
/ = 1, and describes a true physical radiative mode (for 
the electromagnetic field). Eq. (p6|), however, does not 
describe any physical mode, as it is nothing but the evo- 
lution equation for the violation of the momentum con- 
straint, which propagates hyperbolically on its own (cf. 
Eq. (68) of Ref. [^ and the discussion therein). If we as- 
sume the momentum constraint to be satisfied on the ini- 
tial spacelike hypersurface (namely, if we take both Z} 
and its time derivative to vanish on the initial slice), then 
Z[ will vanish forever, and consequently Eq. (pO) car- 
ries no physical information. Consequently, both H^ 
and H2 are proportional to Zj , and Eqs. ( |2^ ) and 
( p5| ) are diffeomorphic and describe the same physical 
radiative mode of the electromagnetic field. The reason 
for the appearance of the equation for the evolution of 
the violation of the constraint equation is clear. The 
field equations wc use are unconstrained, and therefore 
one should not be surprised to find non-trivial equations 
which describe the unphysical sector of the unconstrained 
theory. The field equations we used are the 4>r and the 
(j)9 components of the Einstein equations. In order to 
have a constrained theory, one should include also the 
to and t(j) components, which will insure that unphysical 
modes will not appear. 

Hence, we conclude that in this case there is just one 
radiative mode, as expected from the non-radiative char- 
acter of the dipole gravitational field. We add here that 
in the dipole mode of the polar perturbations (see Paper 
I) all the physical observables (i.e., the metric coefficients 
and the components of the Maxwell tensor) are obtained 
from just one differential equation, even though there 
appear to be two such equations, again as should be ex- 
pected from the nature of the dipole modes. As in the 
axial case we study here, again the unphysical mode is 
an autonomous evolution of the violation of the momen- 
tum constraint. In the axial case one should preclude the 



other differential equation for the evolution of Z} , and 
thus find that the functions H\~' and H2 are diffeo- 
morphic (they are distinguished by just a multiplicative 
factor) . 



B. Decoupling the equation for the general case 

We now turn to the decoupling of the equations in 
the general case (i.e., for any I). We shall follow the 
decoupling procedure of Paper I. We seek functions Z^ 
and Z2 given by the ansatz 

h['^ = aZ[^^ + pz'f^ 



(28) 



and 



H. 



(-) 



lZ{ 



{-) 



5Z. 



2 ' 



(29) 



where a, /3, 7, and 5 are parameters yet to be fixed. Sub- 
stitution into Eqs. (E2h and (|23|) yields 



aK^z[-^ + PK^Z^-^ 



and 



A 



2^(-) 



2g,7 + Z(/ + l)ra + 4-^a 
r 

r 



-A-) 



-A-) 



(30) 



7A-Zi 



SA'Z. 



2y(-) 



[2Q,il-l){l + 2)a 



,Ql 



(-) 



+ 1(1 + l)n - 6M7 + 4-^7]Z| 
r 

+ ^[2Q,{1 - 1)(/ + 2)/3 + l{l + l)r5 

-6MS + 4:-^S]z\~>. 
r 



(31) 



We now multiply Eq. (|30|) by 7 and Eq. (|31|) by a, 
and subtract the two equations. To have the equations 
decoupled we require that 

20*7^ + 6Ma7 - 2Q^{1 - 1){1 + 2)o? = 0. 

This is a quadratic equation in 7, say, whose solution can 
be determined uniquely, once we recall that 7 vanishes 
for ^ =: 1. Hence, we find that 

1 a 



2Q, 



3M - y/9A'P + AQlil - 1){1 + 2) 



1 a 
2Q 



-92, 



(32) 



where 92 = 3M - y/9AP + iQl{l - 1){1 + 2). 



Next we multiply Eq. Q by S and Eq. Q by /3, and 
subtract the equations. To have the equations decoupled 
we now require that 



K-) 



2QJ^ + 6MpS - 2Q41 -l)il + 2)ff 



0. 



Solving this equation for (5, and requiring that 5 will not 
vanish for I — 1. we find that 






or, 



3Af + yJmP + AQl{l-l)(l + 2) 






(33) 



where gi = 3M + y/9M^ + AQKl - 1){1 + 2). The pa- 
rameters a and /3 remain to be fixed arbitrarily, although 
one could let them equal their dipole-case counterparts. 
Substitution of the expressions for the decoupling pa- 
rameters and the ansatz (ESh and (29) in Eqs. ( p2| ) and 
( p3| ) we find that the two decoupled wave equations as- 
sume the form 



A-) _ 



A 



^'^r - ^ 



92 



l)r 



92 -qi 



{2qi - 6M) 



Z- 



(-) 



(34) 



and 



A^zi^^ 



A 



l{l + l)r + A 

qi 



Ql 



qi -92 



(2(72 - 6M) 



^(-) 



(35) 



These equations can be put in a more compact and sym- 
metrical form as 



K^Z^-^ 



v'--^z^-^ 



where the effective potential is given by 



.(-) 



A 



vr' = =T 



1(1 



l)r + 4-^ 
r 



9j 



(36) 



(37) 



and i,j = 1,2, i ^ j. This wave equation is similar 
to the equation given in Ref. [pi. However, in our case 
this equation holds for all the multipole modes, including 
the dipole, whereas in Ref. [^ the dipole modes are ex- 
cluded, because of the inapplicability of the perturbative 
formalism to the dipole modes. 



VI. THE COMPLETION OF THE SOLUTION 

We can now complete the solution be giving the metric 
coefficients and the components of the Maxwell tensor in 



terms of the functions Z^ ' (for dipole modes just Zj ). 
Knowing the functions Z> ' the functions 5| ' can be 
calculated using Eqs. (pS) and (Ea). Then, the functions 
Q and B can be computed by Eqs. (|2^) and (|2^). Using 
Eq. (ph and the (linearized) Maxwell equations (H), m), 
and (0) all the perturbed components of the Maxwell 
tensor are calculable. From Eq. (|2l| ) the function Q is 
known. Then, due to Eq. (y), we know the combination 
92,3 — 93,2 of the perturbations of the metric coefhcients. 
Using Eqs. { {Wj ) and ( |ll| ) we can now obtain explicit 
expressions for the metric coefficients, and thus complete 
the solution. 

We now have a perturbative formalism which is valid 
for the treatment of all the multipole modes, including 
the dipole, and which generalizes the formalism given 
in Ref. (5|. (The formalism for the polar modes is 
given in Paper I.) In a similar way to Ref. M, one can 
now show the relations between the polar and the ax- 
ial modes. (This is done most easily via the Newman- 
Penrose formalism — see Ref. Ipl.) 
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